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SECTION  1 
INTRODUCTION 


1.1  BACKGROUND  INFORMATION 

The  projection  moire  method  has  been  used  successfully  to 
determine  the  out-of-plane  motions  of  large  surfaces  experiencing 
a  dynamic  response  to  an  impulsive  load  (References  1  and  2). 

This  method  will  perform  equally  well  in  situations  where  the 
response  of  the  test  specimen  is  static.  Use  of  this  moire  method 
permits  the  motion  of  the  test  specimen  to  be  determined  without 
physical  contact  between  the  test  specimen  and  the  device.  If, 
then,  surface  strains  can  be  computed  from  the  out-of-plane  moire 
deformations  of  a  test  object,  the  projection  moire  method  would 
become  an  economical,  noncontact  experimental  tool  in  the  stress 
analysis  of  complicated  structural  shapes. 

Determination  of  the  surface  strains  from  the  surface 
deflection  data  presents  some  problems  because  of  error  magnifi¬ 
cation  that  results  from  using  inexact  experimental  data  to 
compute  the  second  derivatives.  Various  numerical  curve  fitting 
methods  have  been  used  to  obtain  the  higher  order  derivatives. 
Among  these,  the  three  methods  used  most  successfully  are  spline 
functions  (References  3  to  7) ,  beam  functions  (Reference  8) ,  and 
least  square  regression  analysis  (Reference  9) . 

1.2  PROGRAM  OBJECTIVE  AND  SUMMARY 

The  primary  objective  of  this  program  was  to  determine 
surface  strains  on  flat  plate  specimens  experiencing  static  loads 
from  out-of-plane  deformations  obtained  with  use  of  the  projec¬ 
tion  moird  method.  To  accomplish  this,  a  four-point  bending  test 
specimen,  in  conjunction  with  a  four-point  bending  load  device, 
was  selected  as  a  target  since  a  fairly  uniform  strain  field 
exists  in  the  midsection  of  the  loaded  specimen.  A  conventional 
method  of  strain  measurement  was  chosen  to  provide  data  for 
comparison  with  the  calculated  moir£  strains.  For  this  purpose, 


1 


a  number  of  electrical  resistance  strain  gages  were  installed  on 
the  test  specimen  and  strain  measurements  made  for  each  load 
condition. 

Theoretical  strains  were  also  determined  to  aid  in  the 
evaluation.  The  small-deflection  beam  theory  was  the  only 
analytical  method  which  permitted  deflections  and  strains  in  the 
four-point  bending  specimen  to  be  predicted  directly.  Therefore, 
a  large-deflection  beam  theory  which  integrates  the  Bernoulli- 
Euler  beam  equations  by  use  of  elliptic  integrals  (References  10 
and  11)  and  a  small-deflection  plate  theory  which  integrates  the 
Lagrange  plate  equation  by  use  of  a  series  of  trigonometric 
functions  (Reference  12)  were  developed  for  use  with  the 
four-point  bending  specimen.  The  moire  deflections  and  strains 
were  compared  with  the  analytically  predicted  deflections  and 
strains,  with  the  experimentally  determined  strains  from  the 
electrical  gages,  and  with  deflections  from  dial  gages. 
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SECTION  2 

EXPERIMENTAL  TECHNIQUES 

A  brief  description  of  the  four-point  bending  test  fixture 
test  specimen,  moire  device,  and  strain  gage  instrumentation  is 
given  in  the  following  paragraphs.  The  tests  performed  during 
this  program  are  summarized  later  in  this  section. 

2.1  TEST  SPECIMEN 

The  test  specimen,  shown  in  Figure  1,  was  a  26-inch  long, 
6-inch  wide,  0.25-inch  thick  flat  plate  made  of  2024-T351  alumi¬ 
num  alloy.  One  side  of  the  plate  was  prepared  for  recording 
moire  fringes  by  the  application  of  a  reflective  coating  and 
installation  of  fiducial  marks  as  illustrated  in  Figure  2. 

2.2  FOUR-POINT  BENDING  LOAD  DEVICE 

A  four-point  bending  load  device  generally  consists  of  two 
knife-edge  or  cylindrical  supports  on  which  a  beam  can  be  simply 
supported  at  its  ends  and  two  loading  points  equidistant  from 
the  end  supports.  When  loaded  this  way,  a  constant  bending  moment 
is  generated  in  the  beam  between  the  loading  points. 

The  four-point  bending  load  device  designed  and  fabricated 
for  this  investigation  is  shown  in  Figure  2.  This  fixture  has 
several  unique  features.  The  beam  is  supported  by  the  inner  load¬ 
ing  points  and  the  loads  are  applied  to  the  ends  of  the  beam  by 
means  of  finely  threaded  power  screws.  The  loading  points  which 
come  in  contact  with  the  test  specimen  are  four  ground  cylindrical 
rods,  3/4  inches  in  diameter  and  6  inches  long.  The  loading  points 
in  the  ends  are  guided  so  that  the  applied  loads  are  parallel  to 
each  other  and  normal  to  the  surface  of  the  test  specimen.  Load 
cells  incorporated  in  the  inner  loading  points  were  calibrated 
prior  to  their  use.  Loads,  i.e.,  analog  voltages,  were  read  on 
separate  digital  voltmeters  to  an  accuracy  of  +.01  volts  or 
+0.5  pounds.  Loads  were  accurately  balanced,  controlled,  and 


6.0  in 
U5.24  cm) 


MATERIAL:  ALUMINIUM  ALLOY 
2024-T351 


Figure  1.  Four- Point  Bending  Specimen. 


Figure  2.  Four- Point  Bending  Load  Device. 


applied  by  turning  the  power  screws  and  monitoring  the  loads  on 
the  load  cells. 

2 . 3  INSTRUMENTATION 

The  moire  device  used  in  this  investigation  is  described 
in  Reference  1.  For  reader  convenience,  the  basic  experimental 
setup  is  shown  in  Figures  3  and  4. 

Four  Micro-Measurements  EA-13-250BF-350  strain  gages  were 

installed  in  the  center  and  back-to-back  on  the  front  and  rear 

surfaces  of  the  test  specimens  in  the  x  and  y  directions  so  that 

they  measured  the  e  and  e  strains  (Figure  5) .  These  strain 

x  y 

gage  circuits  were  balanced  on  a  10-channel  Vishay  bridge  balance 
and  were  read  on  a  Vishay  strain  meter. 

A  dial  indicator  gage  whose  sensitivity  is  0.001  inch  was 
used  to  measure  the  deflections  at  the  center  of  the  plate. 

2.4  -  MATRIX  FOR  EXPERIMENTAL  INVESTIGATIONS 

The  experimental  effort  conducted  in  support  of  the  program 
objective  is  summarized  in  Table  1. 
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Figure  4.  Schematic  View  of  Moir£  Device  Illustrating 
Components  and  their  Relationship  to  Target 
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-SECTION  3 

ANALYTICAL  TECHNIQUES 


The  test  specimen  was  subjected  to  a  range  of  loads  which 
would  induce  both  small  and  large  deflections.  Hence,  it  was 
necessary  to  predict  the  static  response  of  the  test  specimen  by 
both  small  and  large  deflection  theories.  A  brief  outline  of 
the  theories  used  and  the  solutions  obtained  are  given  in  the 
following  paragraphs. 


3.1  SMALL-DEFLECTION  PLATE  THEORY 

The  Lagrange  plate  equation  is  given  by  (Reference  12) 

4  4  4 

3Wj_^  3  w  ^  3  w 

— j  +  2  ~2 — j  +  — j 

3x  3x  3y  3y 

The  simple  support  boundary  conditions  are 


a 

D 


w  *  0 , 


32  w 


0  for  x  *  0  and  x 


3x 


(1) 


(2) 


The  free- edge  boundary  conditions  are 


[  3y^  3x 


0  for  y  ■  +b/2 


(3) 


aiy*  «-*>  4U 

l  3Y  3y  J 


»  0  for  y  ■  +b/2 


The  deflected  surface  is  assumed  in  the  form  of  a  series 
(Figure  5) 


(4) 


w 


W1  +  W2 


(5) 


10 


where 


\  I  l  h  sinT 

A  ir^bD  111*1,3,5  j=l  m 


.  miTX 

sm - 

a 


represents  the  deflection  curve  of  a  four-point  bending  strip. 


4  —  m=l ,3,5 
it  bD 


Am  ==Sh^  * 


inh™]  sinS^ 


A  and  B  are  determined  from  the  boundary  conditions  (2) ,  (3) , 
mm 

and  (4)  and  are  represented  by 


(,1-v)  cosha  -(l+v)  sinha 
m  m 


(\i^+2v-3)  sinh  a  cosh  a  +  {l-\>)  a  1*1 

Tfl  TTI  TTI  » 


mire  . 


- (1-v)  sinh  a 

m 

— 5 - - - 

(v  +2v-3)  sinh  a  cosh  a  +(l-v) 

m  m 


2 

l  sin 
j=l 


nurc  . 


This  series  solution  converges  very  rapidly.  The  first 
three  terms  are  sufficient  for  deflections  and  the  first  six  terms 
should  be  considered  for  stresses.  The  details  of  the  analytical 
expressions  are  given  in  Appendix  A. 

3.2  LARGE-DEFLECTION  BEAM  THEORY 

By  symmetry  the  beam  is  considered  to  be  fixed  at  A  and 
only  one  half  of  the  beam  is  shown  (Figure  6)  .  The  Bemoulli- 
Euler  equations  for  the  beam  in  the  regions  BC  and  AB  are 
given  by 


p 

- - 6. 


(10) 


D~  *  P(x  -x)  +  P  cotij/  (y  -y)  for  x  >  x  >  x_ 
as  C  O  C  C  a. 

D§|  =  P(xc-x)  +  P  cot*o(yc-y) 

-  P (Xg-x) — P  COtX  (y  -y) 

for  x_  ~  x  >  0  (11) 

The  boundary  conditions  are 


g*]-° 

•“  J  ip=4>2 

(12) 

‘*>s-0  “  0 

(13) 

to  be  continuous 

Lds  J«-*1 

(14) 

• 

These  differential  equations  (10)  and  (11)  are  integrated 
and  the  constants  of  integration  are  determined  from  the 

boundary  conditions  (12),  (13),  and  (14).  The  equations 
L^,  L2,  Xg,  yfi,  xc,  and  yc  are  expressed  in  terms  of 
elliptic  integrals  of  the  first  and  second  kind  and  are 
given  below.  The  details  of  the  derivation  are  given 
in  Appendix  B. 

for 

L.  ■  —  K  (p,  ,N) 

A  k 

K1 

(15) 

L2  -  (K (p)  -  k (p,5) ] 

2 

^Pi 

(16) 

yB  -  jq=-  11-  COSN] 

yc  -  yB  +  J-{-sin\|>o  [K(p,m)  -  K(p) 

(17) 

+  2E(p)  -  2E  (p,m)  ]  +  2p  cosi|>o  cosm> 

(18) 

13 


(1?) 


CB  *  jq-  [3K(p1?N)  -  2E(Pl,N)] 


where 


xc  =  Xg  +  JJ—  {cosiJ>0  [K (p,m)  -  K(p) 

+  2E(p)  -  2E(p,m)]  +  2p  siniJ>o  cosm} 


(20) 


*/2 


K(m) 


p  sin  <p 


n/2 


E  (m)  = 


A-p2  sin2$  d$ 


are  complete  elliptic  integrals  of  first  and  second  kind,'  and 


K(p,m)  = 

rm 

/C? 

4 

c 

E(p,m)  » 

3_ 

rm  _ _ _ 

/l-p2 

✓ 

o 

-  .“7* 

-p  sm  $ 


•p  sin  <p  d$ 


are  elliptic  integrals  of  first  and  second  kind, 
p  *  sinir/4 


pi a  757 +  ~2  0 

ki“  /is 


cos(w/2  +  \p 


r*:1]] 


(coti(<o  -  cotXQ)j 


(21) 

(22) 

(23) 


14 


(24) 


m 

N 


D  smf, 


sin-^ 


sin  ^ 


sin  1  o 


P 

sin^/2 


?! 


(25) 

(26) 


A  numerical  solution  of  these  equations  for  a  four-point 
bending  test  specimen  is  obtained  by  applying  the  condition  that 
the  supports  are  a  known  distance  apart  and  the  length  of  the 
beam  can  change  during  loading  (References  10  and  11) .  This 
condition  is  imposed  by  the  use  of  Equations  19  and  20.  A  value 
for  \l>2  is  assumed  and  Equation  20  is  solved  for  m.  The  value 
of  ^  is  then  obtained  from  Equation  25  and  the  value  of  px 
and  N  are  then  computed  from  Equations  23,  22,  and  26  and 
substituted  into  Equation  19.  This  numerical  iteration  is 
continued  and  a  set  of  values  for  4^  and  i|»2  are  finally  obtained 
after  a  few  trials  to  simultaneously  satisfy  Equations  19  and 
20. 


15 


" ,r'  ■  -  -  -  -  * 


Y  infiiia'Tt  -il  ntn 


SECTION  4 

EXPERIMENTAL  AND  ANALYTICAL  RESULTS 

A  brief  description  of  the  experimental  procedure  and  the 
moire,  the  strain  gage,  and  the  dial  gage  results  are  given  in 
the  following  paragraphs.  The  analytical  results  are  presented 
and  compared  with  experimental  results. 

4.1  EXPERIMENTAL  PROCEDURE 

The  total  number  of  tests  were  organized  into  three 
separate  categories,  e.g..  Test  Series  I,  II,  and  III,  as  shown 
in  Table  1.  In  Test  Series  I,  the  moire  device  was  held 
horizontally  thereby  projecting  vertical  Ronchi  rulings.  A 
preload  of  five  pounds  was  applied  to  the  specimen  to  firmly 
position  the  specimen  against  the  supports.  The  strain  gages 
were  zeroed  and  the  reference  moire  fringes  were  recorded.  Loads 
were,  applied  to  the  specimen  in  increments  of  50  pounds  to  a 
maximum  load  of  305  pounds.  At  each  increment  of  load,  two 
photographs  of  the  moire  fringe  pattern  were  taken,  and  the 

meter  readings  were  scanned  manually  and  recorded. 

•  • 

In  Test  Series  II,  the  moire  device  was  rotated  90°  so 
that  horizontal  Ronchi  rulings  were  projected  onto  the  test 
specimen.  The  moire  fringes  and  strains  from  strain  gages  were 
recorded  by  the  procedures  outlined  for  Test  Series  I. 

In  Test  Series  III,  a  dial  gage  was  positioned  to  measure 
deflection  at  the  center  of  the  specimen.  A  preload  of  five 
pounds  was  applied  to  the  specimen  to  firmly  position  the  speci¬ 
men  against  the  supports,  and  the  initial  reading  of  the  dial 
gage  was  recorded.  Loads  were  then  applied  to  the  specimen  in 
increments  of  50  pounds  to  a  maximum  load  of  305  pounds  and  then 
decreased  by  50  pound  increments  to  the  initial  five-pound 
preload.  Dial  gage  readings  were  noted  at  each  increment  or 
decrement  of  applied  load.  The  deflection  at  each  load  was 
obtained  as  the  difference  between  the  dial  gage  reading  at  that 
load  and  the  preload. 
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4 . 2  MOIRE  RESULTS 


The  test  specimen  was  considered  to  be  a  plate  with 
2  2 

curvatures  — ^  and  — £  in  the  x  and  y  directions*  The  moire 
Sx'*  3y 

device  was  used  to  determine  these  curvatures .  The  moire  device 

3^w 

was  held  in  the  horizontal  position  to  determine  — y  and  in  the 

32w  3x 

vertical  position  to  determine  — y.  Moire  fringe  patterns  for 

these  cases  are  shown  in  Figures  7  and  8.  Deflected  shapes  of  the 
loaded  plates  were  computed  using  measurements  taken  from  the 
photographs  and  a  program  called  STRAIN.  Deflection  data  were 
fitted  to  quadratic  and  cubic  polynomials  using  a  least-square 
curve  fitting  routine  called  LSQ.  Slopes  and  curvatures  of  the 
deflected  shapes  were  determined  using  a  program  called  PLATE. 

The  moire  strains  are  given  from  the  relations 


e  a 

X 


In  the  case  of  large  deflections,  the  curvatures  are  replaced 
by  the  expression 


i  +  tel2] 

{  3xJ 


[■  *  (sff 
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(a)  5  POUNDS  LOAD 


<b)  55  POUNDS  LOAD 


(c)  105  POUNDS  LOAD 


(d)  155  POUNDS  LOAD 


<f)  255  POUNDS  LOAD 


<q)  305  POUNDS  IOAD 


(b)  55  POUNDS  LOAD 


(c)  105  POUNDS  LOAD 


Figure  8, 


Moire  Fringes  to  Dote  mono 


Curvature , 


20 


(e)  205  POUNDS  LOAD 


(g)  305  POUNDS  LOAD 


Figure1  8.  (Concluded) 


2  l 


3^w 

The  — r  cvirvature  of  the  four-point  bending  test  specimen 
3y  32w 

is  a  small  quantity  compared  to  — *•,  and  the  very  small  number  of 

3x 


moire  fringes  are  insufficient  to  accurately  determine  this 

Quantity.  Therefore,  according  to  References  13,  14,  and  15, 

*2  2 
3  w  3  w 

— j  curvature  is  assumed  to  be  -v — *•  and  the  strain  e  is  given  by 
3y  3x  y 


e 

y 


-V 


e 


x 


(31)  • 


The 


32w 


curvatures  are  also  determined  directly  from  Figure  8 


3y 


for  a  few  load  cases. 


The  moire  deflections  and  strains  are  presented  for  loads 
of  50  to  300  pounds  in  Tables  2  and  3,  and  in  Figures  9  and  10. 

The  moire  deflections  computed  by  program  STRAIN  were  fit  to  a 
.  second  order  polynomial  using  the  least  squares  routine  PLSCF 
and  a  second  order  Legendre  polynomial,  and  the  experimental 
points  and  the  fitted  curves  were  plotted  using  plot  routine 
PLOTS,  Figures  11  and  12. 

The  scatter  of  moire  experimental  data  is  large  at  small 
loads  and  the  experimental  data  points  cannot  be  fit  by  means  of 
a  polynomial.  It  should  be  remembered  that  the  plate  specimen  has 
initial  imperfections  which  are  of  the  same  order  of  magnitude 
as  the  deflections  corresponding  to  these  smaller  loads.  This 
could  be  the  cause  for  the  large  amount  of  scatter.  At  higher 
loads,  the  moird  experimental  points  can  easily  be  fitted  by 
means  of  second  or  third  degree  polynomials  and  the  experimental 
scatter  is  greatly  reduced. 


4.3  DEFLECTIONS  FROM  DIAL  GAGES 

The  deflections  measured  by  dial  gages  are  given  in 
Table  4  and  shown  previously  in  Figure  9. 
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TABLE  2 

MOIRE  DEFLECTIONS  AND  STRAINS 
AT  THE  BEAM  CENTER 


Load 

(lbs) 

Deflections 

W 

(in) 

Strain 

(U) 

ex 

ey 

50 

0.045 

385.3* 

-128.44 

(363.64)+ 

(-121.22) 

100 

0.094 

765.84 

-255.28 

(766.34) 

(-255.46) 

150 

0.144 

1185.90 

-394.9 

(1186.66) 

(-394.9) 

200 

0.193 

1599.18 

-533.06 

(1602.68) 

(-534.22) 

250 

0.251 

1987.76 

-662.58 

(2005.52) 

(-668.52) 

300 

0.300 

2347.38 

-782.46 

(2398.84) 

(-818.30) 

♦Quadratic  polynomial  curve-fit 
tCubic  polynomial  curve-fit 


Wt 


r 


Ween 


.Wp 


TABLE  3 


MOIRE  DEFLECTIONS  AND  STRAINS 
THREE  INCHES  BELOW  CENTER 


Load 

(lbs) 

Deflection 

w 

(in) 

Strain  (u) 

E 

X 

e 

y 

50 

0.061 

465.38* 

-155.12 

(455.3)+ 

(-151.76) 

100 

0.099 

719.76 

-239.52 

(718.56) 

(-239.52) 

150 

0.147 

1166.5 

-388.82 

(1161.9) 

(-387.3) 

200 

0.193 

1542.78 

-514.24 

(1537.32) 

(-512.44) 

250 

0.252 

1981.46 

-660.5 

(2005.54) 

(-668.52) 

300 

.  0.304 

2347.38 

-784.32 

(2406.96) 

(-802.32) 

•Quadratic  polynomial  curve- fit 
tCubic  polynomial  curve- fit 


^cen 


P 


P 
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G  4, 


mm  2.0  |  —  i  mm  2.0 


OTE:  Limit  of  accuracy  on  moire* data  is  0.25mm 

Figure  11.  Regression  Analysis  Using  a  Second-Order  Polynomial  Relating 
Deflections  and  Span  Positions. 


0>  ujui 


Figure  11.  (Continued) 


mm  4  0  - - -  n»w  4.0 

3.5  V  35 


Figure  12.  Regression  Analysis  Using  a  Second-Order  Legendre  Polynomial 
Relating  Deflections  and  Span  Positions. 


Figure  12.  (Concluded) 


TABLE  4 

DIAL  GAGE  DEFLECTIONS  AT  CENTER  OF  BEAM 


Load 

(lbs) 

Deflections 

(in) 

50 

0.052 

100 

0.1055 

150 

0.1590 

200 

0.212 

250 

0.266 

300 

0.322 

P  P 


W  ■  WCEN-WP 


4.4  STRAINS  PROM  ELECTRICAL  RESISTANCE  GAGES 

The  strains .  recorded  f  roih  strain  gages  are  presented  in 
Table  5  and  shown  previously  in  Figure  10. 

4 . 5  ANALYTICAL  RESULTS 

The  deflections  and  strains  predicted  from  the  small- 
deflection  beam  theory  are  presented  in  Table  6 ;  from  the  small- 
deflection  plate  theory  in  Table  7;  and  from  large-deflection 
beam  theory  in  Table  8. 


34 


w 

«  u 

w  ^ 

J  tn 

(0  M 

<  cn 
e  64 
06 


CN 

in 

m 

in 

>i 

• 

• 

• 

H 

W 

r-4 

o 

VO 

35 

1 

CN 

CN 

VO 

r* 

VO 

4p 

>. 

H 

r4 

CN 

cn 

n* 

in 

VO 

u 

>- 

U 

4H 

4H 

4H 

■H 

4H 

4 

■1 

i  a. 

<u  — 

O'  <N 

10  >i 

O  u 


m 

i  a 
<0  — 
o>  «-i 
•(0  >. 
O  w 


CN 

CO 

o 

•*r 

r> 

m 

r* 

ON 

ON 

00 

CN 

m 

N* 

in 

VO 

m  m  m  «w  rH 

n  n  to  *o 

i  i  i  i  i 


PM 

II  X 

<  w 

X  + 

W  r-M 

X 

w 


CN 

1 

VO 

• 

O' 

00 

CN 

VO 

in 

1 

PM 

in 

00 

r«4 

1 

I 

1 

I 

4? 


PM  — . 

i  a 
<D  ~ 
O'  PM 
IQ  X 
u  u 


N» 

CN 

VO 

in 

CN 

in 

|  CN 

N* 

t*> 

fH 

© 

© 

00 

CN 

r* 

m 

i 

1 

H 

•H 

CN 

PM 

i 

i 

•i 

i 

1 

o 

cn 

o 

fH 

© 

in 

CN 

CN 

CN 

© 

VO 

N* 

00 

CN 

© 

ON 

cn 

+ 

+ 

PM 

+ 

+ 

+ 

+ 

o  ©  o  o  o  o 

in 

o 

in 

o 

in 

O 

rH 

i—4 

CN 

CN 

cn 

- 

TABLE  6 

DEFLECTIONS  AND  STRAINS  FROM  SMALL- DEFLECTION  BEAM  THEORY 


Load 

(lb) 

WLOAD 

(in) 

WCEN 

(in) 

W 

(in) 

£ 

X 

Strain 

u 

50 

0.1638 

.2253 

.0615 

427 

-142.33 

100 

6.3277 

.4506 

.1229 

853 

-284.33 

150 

0.4915 

.  6759 

.1844 

1280 

-426.67 

200 

0.6554 

.9012 

.2458 

1707 

-569 

250 

0.8192 

1.1265 

.3073 

2133 

-711 

300 

0.9830 

1.3518 

.3688 

2560 

-853.33 

TABLE  7 

DEFLECTIONS  AND  STRAINS  FROM  SMALL-DEFLECTION  PLATE  THEORY 


Load 

W 

WLoad 

inch 

w 

CEN 

inch 

W 

inch 

ex 

Strain 

y 

ey 

Strain 

u 

50 

0.1522 

0.2064 

0.0542 

346 

109 

100 

0.3045 

0.4129 

0.1083 

693 

217 

150 

0.4567 

0.6193 

0.1626 

1040 

326 

200 

0.6090 

0.8257 

0.2168 

1386 

435 

250 

0.7612 

1.0322 

0.2710 

1733 

543 

300 

0.9135 

1.2386 

0.3252 

2080 

652 

TABLE  3 

DEFLECTIONS  AND  STRAINS  BY  LARGE -DEFLECTION  BEAM  THEORY 


Load 

W 

W,.  , 

Load 

inch 

W 

CEN 

inch 

W 

inch 

e 

X 

Strain 

u 

e 

y 

Strain 

u 

100 

0.3330 

0.4577 

0.1247 

855 

-285 

200 

0.6760 

0.9308 

0.2547 

1724 

-575 

300 

1.0351 

1.4182 

0.3830 

2621 

-874 

w*  Wcen-Wp 


SECTION  5 

ANALYSIS  OF  RESULTS 

Deflections  are  compared  in  Figure  9  and  in  Table  9,  and 
strains  are  compared  in  Figure  10  and  in  Tables  10,  11,  and  12. 

The  deflections  from  dial  gages  and  strains  from  resistance  gages 
are  considered  as  standards  with  which  the  deflections  and  strains 
by  other  methods  are  compared  and  evaluated. 

In  comparison  with  the  dial  gage  deflections,  the  moire 
deflections  are  smaller  by  7  to  13  percent?  deflections  predicted 
by  small-deflection  plate  theory  are  higher  by  1  to  4  percent? 
the  large-deflection  beam  theory  deflections  are  higher  by 
18  to  20  percent?  and  small-deflection  beam  theory  deflections 
are  higher  by  15  to  18  percent.  Thus  the  small-deflection  plate 
theory  predicts  deflections  in  close  agreement  to  the  dial  gage 
deflections.  The  deflections  determined  by  the  moird  fringe 
technique  are  more  accurate  at  higher  loads.  This  is  due  to 
initial  imperfections  of  the  test  specimen  which  will  be  of  the 
same  order  of  magnitude  as  the  deflections  at  small  loads.  This 
causes  considerable  scatter  of  the  experimental  data  points  as 
seen  in  Figures  11  and  12.  At  higher  loads,  when  the  deflections 
are  large,  the  scatter  in  experimental  data  points  is  reduced  and 
moir4  deflections  agree  within  6  percent  with  dial  gage  deflec¬ 
tions.  Another  cause  for  the  large  discrepancy  noticed  between 
moire  and  dial  gage  deflections  was  due  to  the  fact  that  the 
initial  reference  surface  was  not  perfectly  flat,  as  most  machined 
flat  plates  have  some  initial  imperfections.  The  test  specimens 
were  prepared  from  rolled  plate  and  developed  imperfections  during 
machining  which  could  not  be  easily  removed. 

In  comparison  with  resistance  strain  gage  strains,  the 
moir6  ex  strains  computed  by  second  and  third  order  polynomials 
are  smaller  by  0  to  8  percent  (Table  10  and  Figure  10) ,  whereas 
those  computed  by  Legendre  polynomials  are  smaller  by  6  to  25 
percent  (Table  11) ?  the  classical  beam  theory  tx  strains  are 
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higher  by  1.5  to  3  percent;  the  large- deflection  beam  theory  ex 
strains  are  higher  by  3  to  6  percent;  the  £x  strains  from  classical 
plate  theory  are  smaller  by  16  to  18  percent.  Thus  the  moire 
strains  computed  by  second  and  third  order  polynomials  are  in 
close  agreement  with  resistance  gage  strains,  especially  at 
higher  loads,  and  can  be  predicted  by  small-deflection  beam  theory 
over  the  range  of  loads  considered  in  this  investigation. 

The  moire  e  strains  were  determined  by  two  methods.  In 
the  first  method,  since  the  curvature  in  the  y-direction  is 
(References  13,  14,  and  15) : 


the  ey  strain  can  be  approximated  by 


In  the  second  method, 


curvature  is  determined  directly  from 


the  Moire  fringe  photographs  of  Figure  8.  The  ey  strains  obtained 
by  the  first  method  are  compared  in  Tables  11  and  12  and  in 
Figure  10.  In  comparison  with  resistance  gage  strains,  the 
moire  strains  obtained  by  second  and  third  order  polynomials  are 
higher  by  7  to  34  percent;  whereas  those  obtained  by  a  second-order 
Legendre  polynomial  are  higher  by  6  to  20  percent;  the  small- 
deflection  beam  theory  strains  are  higher  by  15  to  32  percent; 
the  large” deflection  beam  theory  ey  strains  are  higher  by  17  to 

35  percent.  Classical  plate  theory  predicts  tensile  e~ 
strains  and  cannot  be  compared  with  resistance  gage  strains. 

The  solution  strategy  used  in  the  classical  plate  theory  is  similar 
to  the  one  used  by  Timoshenko  for  many  of  the  plate  bending  problems. 
But  this  strategy  may  not  give  true  solutions  due  to  the  type  of 
loading,  boundary  conditions,  and  possibly  due  to  the  large  aspect 
(a/b)  ratio  of  4  of  the  test  specimen.  The  »oir4  strains 
determined  by  the  second  method  for  loads  of  150,  200,  and  250 
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pounds  are  compared  in  Figure  10  and  these  are  smaller  than 
resistance  gage  strains  by  0  to  8  percent.  Therefore,  none  of 
the  analytical  methods  considered  in  this  investigation  predict 
the  e  strains  accurately.  The  e^,  Moire  strain  determined  by  the 
second  method  is  closer  to  resistance  gage  strains  and  therefore 
is  to  be  preferred  to  the  first  method. 

For  the  loads  considered  in  this  investigation,  curvatures 
determined  by  large  deflection  equations  (Equations  29  and  30) 
did  not  differ  significantly  from  the  curvatures  determined  by 
small  deflections  and  therefore  the  numerical  results  are  not 
presented  for  these. 


SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  conclusions  are  based  on  the  results  of 
this  program. 

1.  The  projection  moire  method  used  in  this  investigation 
determines  both  the  deflections  and  strains  accurately  in  flat 
surfaces  undergoing  large  deformations  due  to  static  loads. 
However,  the  strains  determined  by  regression  analysis  using 
second  and  third  order  polynomials  and  second  order  Legendre 
polynomials  are  accurate  only  for  the  center  portion  of  the 
beam. 

2.  Of  the  analytical  methods,  small-deflection  beam 
theory  and  large-deflection  beam  theory  predicted  beam  deflec¬ 
tions  and  ex  strain  accurately.  Plate  theory  predicted  the 
deflections  very  accurately.  The  strains  were  not  predicted 
accurately  by  any  of  the  methods  investigated. 

The  following  recommendations  are  made.  (1)  Other  numerical 
curve  fitting  methods  such  as  spline  approximations  and  beam 
functions  should  be  tried  to  obtain  accurate  second  derivatives 
from  the  moire  deformations  at  any  position  along  the  span  of 
the  beam.  (2)  The  investigations  reported  herein  should  be 
extended  to  evaluate  the  accuracy  of  the  method  for  curved 
surfaces  under  dynamic  response  conditions. 
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APPENDIX  A 


SMALL  DEFLECTION  PLATE  THEORY  FOR  FOUR- POINT 
BENDING  PLATE  SPECIMEN 


BASIC  EQUATIONS 


The  Lagrange  plate  equation  is  given  by 

3  4w  .  „  3  4w  .  3  4w  a 

— I  +  2 — 2 — ?  +  — T  =  - 
3x  3x  3y  3y*  D 


(A-l) 


where  w  is  the  downward  deflection,  q  is  the  intensity  of  the 

Eh  ^ 

load  and  D  is  given  by  - «—  .  The  boundary  conditions  are 

12 (1-v  ) 


w  *  o  ; 


2  2  1 
3  w  .  3 w  _  0 

- y  *r  V - —  U 

3y2  3x2 


+  (2-v) 


3x  3y 


for  x  «  0,  and  a 


(A-2) 


for  y  *  ±b/2 


(A- 3) 


for  y  *  ±b/2 


(A- 4) 


The  load  is  symmetrical  with  respect  to  the  x-axis 
(Figure  4)  and  therefore  the  deflection  surface  is  symmetrical 
with  respect  to  x-axis.  Therefore  we  can  consider  only  the 
conditions  along  the  side  y  »  +b/2. 


The  deflection  surface  is  assumed  in  the  form 


W  *  Wj_  +  Wj 


(A-5) 


where 


,_,3  »  1  1  m*x 

-  Hs  t  t  K  —r1  •  ““V 


*  bD  m-1,3,5  j«l  ® 


w  -  ss-  r  a  coah^B J^-sinh^  sin^ 
2  -4.  7T  4-  m  a  m  a  a  a 

*  bD  in-1,3,5  1 


The  constants  and  B_  are  given  by 
mm 


-  (l-v)coaha  - (l+v)sinha  « 

.  _  ra  in  .  \  I 

m  (v  *2v-3) ainha  cosha  +(l-v)  a  j«l  * 
ram  in 


Bm  “  nr 


- (1-v) ainha 


(v  +2v- 3) ainha  coaha  +<l-v)  a  j-l 

at  in  n 


2  , 
r 

•  <-  ain — — * 


Where  a  *  21k 
m  a 

The  moment- curvature  relations  are  given  by 


—  w 
M  -  -D 

x  L3x 


My  -  -D 


B?  *  -9 


M _ -  D(l-V) 


3x3y 


The  maximum  stresses  at  the  surface  of  the  plate  are  given  by 


The  maximum  strains  at  the  surfaces  of  the  plate  are  given  by 


£  = 


32w  - 

h 

92w 

+ 

z 

=  + 

3x 

2 

•  z? 

32w  - 

h 

32w 

T 

z 

— I  + 

3y 
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’i? 

(A-16) 

(A-17) 


The  value  of  the  first  three  constants  of  the  trigonometric 
series  are  given  by 


B1 

A3 

B3 

A5 

B5 

+0.03558502 

0.01535823 

0.01342625 

-0.0001234 

-0.00012857 
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APPENDIX  B 


LARGE  DEFLECTION  BEAM  THEORY  FOR  A  FOUR-POINT 
BENDING  BEAM  SPECIMEN 


BASIC  EQUATIONS 

El 

Let  L  =  L.  +  L,  be  the  length  of  cantilever  and  D  =  - *• 

1  1-v 

its  flexural  rigidity.  For  a  point  U  between  B  and  C  (Figure  5) 

the  bending  moment  equals 


Mu-  P(xc-x)  +  P  co.t*o(yc-y)  »  d|| 


(B-l) 


where  is  the  angle  between  the  slope  and  the  horizontal,  and 
s  is  the  arc  length  from  A  to  U  measured  along  the  elastic 
shape.  Differentiating  with  respect  to  s,  we  get 


D^-i  =  “P  “  P  cotifc  ^  =  -P  cos ip-P  cotip  sirup 
•  as  o  as  o 


(B-2 ) 


In  terms  of  new  variables 


u  *  EJ  ;  9  «  *  +  *0 


(B-3) 


the  Equation  2  reduces  to 

+  C*sin0  -  0 

du4 

where  C*«  kj2  and  kj2  -  gjf . 

Further  Equation  4  can  be  written  in  the  form 


(B-4) 


7  H  ’  c*«*e  +  CI 


(B-5) 


For  the  solution  of  the  constant  C^  we  have  the  boundary 
condition 


(B-6) 
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(B-7) 


Therefore  «  0,  and 

d9 

3u 


/  2C*cos6 


This  equation  can  be  integrated  to  give 


/L  ,  rif/2 

ds  =  —  | 

l  2SC*  > 


J1  . 


_ d9 

/  sin2Tr/4-sin29/2 


(B-8) 


To  bring  this  Equation  8  to  Legendre's  standard  form,  let 


P  *  sinir/4; 


p  sin<fr  =  sin9/2 


(B-9) 


In  terras  of  the  new  variable  <J>,  the  Eauation  8  reduces  to 
,  n/2  « 

t  ■  pi  i  f  - •» —  -  [  f*  n — 1/2“+  [K(p)-K(p,m) ] 

2  [l-p2,inV/2  l  tl-pW*]1/2  * 

where  m  ■  Sin-  I—  Sin  — 2~^J 

The  equation  for  the  length  s  to  the  point  u  is  given  by 

SU  *  L1  +  jq  (K (p,n) -K (p,m) } 


(B-10) 


where  n  »  Sin”1 


IP  *  Sin 


(B-ll) 


(B-12) 


The  equation  for  the  horizontal  distance  x  to  the  point  d  is 
given  by 


/•n 


dx 


ds  cosij; 


(B-13) 


m 


Integrating  this  equation,  we  get 
xO  *  *8  +  f  eo,*o  tK<P»»>-Mp,n)*2E(p,n)-2Etp,mH+2p  *in*o (co«i-co«n) }  (3-14) 


The  horizontal  distance  to  the  end  of  the  cantilever  C  can  be 
found  by  substituting  an<2  n=ir/2  in  Equation  14. 


Xg+JJ-  {COS*o 


[K(p,m)-K(p)+2E(p)-2E(p,ra) ]+2p  sin*o  coam} 


(B-15) 


The  equation  for  the  vertical  distance  y  to  the  point  U  is  given 
by 


f%y 


rn 

ds  sin 4) 

. 

m 


(B-16) 


Integrating  this  equation, 


{2p  co»*o<coai»-coan)  -»in*o  [K  (p,n)  -k  (p,n)  +2E  (p,n)-2E(D,m)  ) } 


(B-17) 


The  vertical  distance  to  the  end  of  cantilever  C  is 


yB+jj-  {2p  coapgcoam-ainpo [K (p, m)  -K  (p)-+2E  (p)  -2E  (p ,m) 


1) 


(B-18) 


For  a  point  Q  between  A  and  B,  Figure  3,  the  bending  moment 
equals 


“q  -  P(xc-x)*P  cot*o(yc-y)-P(xB-x)-P  coU0<yB-y>  -  (B-19) 


Differentiating  with  respect  to  s,  we  get 


-p  «**<,  a^p  g  -  sin, 


where 


2 


(cotij>  -cotX  ) 
•  o  o 


p 

ff- 


(B-20) 
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Equation  20  dan  be  written  as 


r  \  2 

lfdtlM  ,  2 

l[dij  “  kl  cos*  +  C2 


(B-21) 


For  the  solution  of  the  constant  C2  we  have  the  boundary 
condition 


mi  =  /2  k2  /cbs  (iK+'J)  ) 


Solving  for  C 2 


C2  =  2[k22  cos  ('l'1+'('0)  -  h12  costal 


(B-22) 


Now  the  Equation  21  becomes 


/2  k^  (cosi|»-cosi|;1)+— cos  (♦1+<PQ  )  j 


(B-22) 


Introducing  a  new  variable  such  that 


{cos  +  cosi^ 


(B-23) 


we  can  integrate  the  equation 


Li  _  i  f*i  <tt 
a*  1ST  7=f^ 

J  Ai  *sin 


(B-24) 


To  bring  this  Equation  24  to  Legendre's  standard  form,  let 


sin4>/2  *  p,  sinQ 


(B-25) 


Then 
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L  ^  ~ 
1  *]. 


N 


de 


/l-p^sin2^ 


k7k(Pi,N) 


-i  i 

where  N  =  sin  [ —  sin— s-] 

Pi  2 

The  horizontal  distance  x  to  a  point  Q  is  given  by 


rxn  /N 

u  dx  *  u  ds  cos^ 


XQ  “  3 


ae 


/l-p  2sin2 


p.'sin'O  o 


2 |Nu  A-Pl2sin28  de  -  i-  (3fk(p1,Nu)-2E(p1,Nu) > 


where  N  »  sin-1  t~  s^nT3 
u  p,  z 


The  horizontal  distance  Xg  to  the  point  Q  is  given  by 
Xg  »  i-  |3[k(p1,N)  -  2E(Pl,N)]} 

.  -1  T  1  *1  1 

where  N  *  sxn  Jjj—  sin-j  J 


The  vertical  distance  y  to  a  point  Q  is  given  by 

ry 

ly  *  |  ds  sin*  *  —i-  [1-N,J 

kl 


rw 


k. 


The  vertical  distance  YQ  to  the  point  B  is  given  by 

2pl 

*  TET  tl-NI 


(B-26) 


(B-27) 


(B-28) 


(B-29) 


(B-30) 


(3-31) 
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